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baiidwidth  to  100  >U;  (hu&  the  ballscrow  was  obie  to  reach  a maximum  accoleraiion  of 
2 g.  The  complexity  in  designing  digital  rillers.  their  long  calculation  time  and  their  lack 
of  robustness  make  them  undesirable  in  a machine  tool  controller.  Therefore,  a sliding 
mode  controller  was  implemented  as  a substitute  for  the  PD  controller  with  digital  filters 
to  acquire  a high  bandwidth.  ITus  robust  controller  permitted  a bandwidth  of  65  Hz  to 
achieve  an  occclcmtion  of  I g. 

A linear  motor  drive  under  PD  control  was  found  to  have  a 70  Hz  bandwidth.  No 
additional  filleni  were  requited  to  achieve  this  bandwidth.  A major  drawback  of  this 
drive  is  its  sensitivity  to  changes  of  mass  and  to  exteraal  disturbances  like  cutting  foires 
and  friction. 

Feedforward  filters  were  applied  on  both  drives  to  decrease  tracking  erroni  during 
high  accelenition  moves.  These  feedforward  fillers  nrduced  the  tracking  error  to  four 
counts  on  the  ballscrew  drive  and  to  two  counts  on  the  linear  motor  drive. 


CHAPTER  1 
INTRODUCTION 

High'.speed  machining  evolved  in  die  paal  decade  aa  an  improvement  lo 
convcnilonaJ  machining.  Advancements  in  microprocessors  and  digital  signal  processors 
(DSPs)  made  it  possible  for  the  computer  numerical  control  (CNC)  to  accurately  control 
feed  drives  at  high  speeds.  The  processor  evolution  coupled  with  a high  demand  on 
delivering  machined  pans  at  lower  cost  and  m Taster  limes  prompted  the  machine  tool 
industry  10  produce  machining  centers  dedicated  only  lo  high-speed  machining. 
Furthermore,  high-speed  machining  made  it  possible  to  replace  other  traditional 
manufxluring  processes.  For  caampic  high-speed  machining  of  aluminum  made  it 
possible  to  produce  monolithic  pans  to  replace  parts  robricaied  from  several  components. 
A machined  part  con  be  more  accurate  and  structurally  siirfcT  than  a riveted  pan.  From 
Zelinski  [i].  the  speed  bralte  of  an  F-15  nghterjel,  made  by  Boeing  Aerospace,  can  now 
be  machined  from  one  piece  of  solid  aluminum.  Traditionally  the  1 1 fool  long.  3 foot 
wide  pan  was  made  from  sheet  meialt  a process  that  lakes  three  months  to  assemble  all  the 
500  components  for  ore  brake.  The  whole  pan  can  now  be  machined  in  3U  hours  on  a 
high-speed  machine.  The  lime  and  cost  of  the  pan  is  uemendousiy  reduced,  but  more 
imponandy,  the  quality  of  the  pan  Is  improved. 

Another  exompie  is  the  high-speed  machining  of  hiird  steel  dies.  High-speed 
machining  ofdies  diminales  some  of  the  steps  from  the  imdiiional  manufacturing  process 
for  dies.  High-speed  machining  makes  it  possible  for  the  dies  lo  be  machined  by  taking 


smaller  loci  steps  Ihus  improving  surface  finish  and  dirainaiing  grinding  and  hand 
polishing. 

Several  aspects  make  a machining  ceniersuilahle  for  high  speed  machining: 
maximum  spindle  speed,  maximum  feed  drive  speed,  lookahead  and  high  dole 
iransmiasion  rales.  Lookahead  is  a CNC  feaiurc  that  analyzes  ihe  commands  beforehand 
10  determine  if  the  machine  is  capable  of  handling  the  commanded  moves.  Smith  [2] 
provided  a more  deiailcd  de.scriplion  of  high-speed  machine  design. 

Although  all  of  Ihe  above  aspecis  are  necessary  for  a high-speed  machine,  this 
work  concenlraies  on  ihe  feed  drive.  Spocincnlly,  ihe  difierem  lypes  of  feed  drives  that 
are  suitable  for  these  machines  and  the  control  methods  ihni  are  implemented  to  make 
them  perform  coordinated  moves  at  high  speeds. 

Feed  drives  on  high-speed  machining  centers  can  exceed  speeds  uf  0.5  m/s.  hut 
that  alone  does  not  make  them  truly  high  speed.  High  acceleration  is  required  in  order  to 
attain  the  maximum  velocity  within  the  commanded  move.  A high-speed  feed  drive  must 
have  high  ucccicraiion  and  high  velocity. 

High  Speed  Macltining 

Most  machined  oircrah  components  like  those  in  Figure  1-1  have  small  pockets  to 
reducelheir  mass  without  diminishing  their  strength.  These  pockets  require  that  Ihe  feed 
drives  of  ihe  machines  reach  (he  required  speed  as  fast  as  possible  in  order  to  produce  cost 
effective  pans.  Feed  drives  on  eonvemional  machines  can  lypically  accelerate  up  to  0.3  g 
which  is  usually  not  sufllcienT  For  high-speed  machining  of  pockets.  Higher  acceleration 
is  requiied  to  lake  full  advantage  of  the  feed  drives  over  short  distances.  An  example  of 
Ihe  need  for  hi^cr  acceleration  is  die  approximation  oF  a helicopter  rotor  grip  in 


Figure  l*t.  Example  of  chin  aJuminum  pans 


Figure  I -2.  Il  y/sa  mode  from  Al  7075-T6  and  was  machined  using  a convenlional  and  a 
high-speed  machine.  The  main  pockei  in  ihc  pan  was  63  mm  by  1 14  mm.  In  Case  A ihc 
conveniionai  machine  commanded  a feedrate  of  25.4  m/min  aiU.2  g over  the  114  mm 
travel.  It  took  the  drive  220  ms  to  acccicraui  to  the  commanded  speed;  then  it  traveled  for 
45  ms  at  constant  speed  before  it  started  its  deceleration.  Most  of  the  move  was  spent 
accelerating  and  decelerating  the  drive.  In  Case  B the  high-speed  machine  was 
commanded  the  some  position  and  velocity  but  a 2 g acceleration.  The  drive  reached  its 
masimum  speed  in  22  ms  and  traveled  at  the  musiroum  speed  for  247  ms  before  .slowing. 
The  total  move  was  performed  in  29 1 ms  while  it  look  485  ms  incase  A.  The  entire  pan 
with  iu  other  features  took  53  minutes  to  be  machined  on  a conventional  machine  and 
look  only  12  minutes  lobe  high-speed  machined.  In  the  amount  of  lime  it  took  to  produce 
one  pan  on  the  conventional  machine,  five  pans  could  be  produced  on  the  high-speed 
machine. 

In  summary,  Il  is  important  to  increase  the  average  velocity  in  order  to  reduce  part 
production  lime.  The  goal  is  to  achieve  the  hipest  acceleraiion  possible. 

Acceleration  and  Bandwidth 

Acceleration  is  limited  by  two  factors:  thrust  capacily  of  the  drive  and  machine 
dynamics.  Motors  have  limited  peak  accelerations  defined  by  ihe  maximum  curreni  that 
can  be  applied  across  the  windings  at  a given  time.  In  general,  this  is  nonhe  limiting 
factor.  Acceleration  is  usually  limited  by  machine  dynamics.  Vibrations  in  the  drive 
preveni  the  machine  from  reaching  high  accelerations.  The  reason  is  that  vibrations  lirail 
the  overall  posicionai  bandwidth  of  the  system.  The  bandwidth  of  the  system  is  an 
indication  of  how  the  system  responds  lo  imnsienis  in  Ihe  lime  domain. 
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Figure  1-3  shows  how  bandwidih  and  acceleration  are  related.  From  Cased  undfi 


in  Figure  1-2.  Ihc  desired  accelerations  look  like  impulse  commands.  The  Fourier 
transforms  ofthe  impulses  show  the  required  bondwidths.  The  acceierations  for 
mxhining  the  rotor  grip  require  two  different  system  bandwidth.  Cased  requires  a very 
low  bandwidth  while  Case  ft  requires  10  times  mote  bondwidths.  It  can  be  .seen  that  high 
acceleration  requires  high  bandwidth.  The  challenge  in  attaining  high  bondwrdth  is  the 
attenuation  of  structural  modes  within  ihe  desired  bandwidth.  Conventional  ballscrew 
drives  generally  have  inodes  less  than  100  Hz.  These  modes  requite  low-pass  rilicnng  to 
reduce  their  e%cts.  These  low-pass  filters  limit  the  overali  bandwidth. 

IVo  approaches  can  be  taken  to  improve  bandwidth:  the  first  is  to  re-design  the 
drive  to  have  all  modes  above  100  Hz  or  to  use  a different  conirol  method  that  attenuates 
the  modes  in  the  control  loop.  The  first  approach  Is  the  most  widely  accepted  in  industry. 
Conventionally  designed  drives  with  rotary  motors  and  leadscrews  are  being  replaced 
with  linear  motors.  Because  linear  moiors  usually  have  modes  at  higher  frequencies  they 
am  suitable  substitutes  for  ballscrews  in  some  cases.  They  have  iimiung  factors  that  are 
different  than  ballscrew  drives  and  cannot  be  switched  without  proper  analysis.  The 
second  approach  is  lo  modify  the  control  system.  Control  methods  can  attenuate  or  reject 
modes  and  lncrea.se  Ihe  overall  bandsvidth.  ThedrRiculty  is  In  implemenling  these  control 
schemes.  Off-the-shelf  controllers  ore  programmed  for  PID  conirol  which  is  a simple 
linear  control  theory  that  is  widely  understood  and  used;  but  this  theory  is  not  capable  of 
dealing  with  structural  modes  within  the  bandwidth.  A more  elaborate  conirol  method  is 
required. 


Figure  1-3.  Acceleralion  and  bojidwidth 


Finally,  moving  feed  drives  al  high  speeds  develops  large  (rocking  errors. 

Tracking  errors  are  (he  difference  be(ween  (he  commanded  and  uciuol  posidons,  as  shown 
in  Figure  1*4.  These  errors  con  cause  imperfecdons  in  (he  machined  pans,  especially 
when  (here  are  muUlpte  aaes  moving  simuUaneously.  Again,  con(rol  (hcoiy  can  reduce 
(hese  errors  and  improve  pan  qualiry, 

Tesi  Beds 

The  Machine  Tool  Research  Cenler(MTRC)  ai  (he  Universlly  of  Florida  has 
designed  (wo  high-speed  milling  machines  (o  be  used  as  rest  beds  for  (he  dcvciupmcni  of 
high-speed  machining.  Figure  I -5  shows  (he  UF-MTRC  HSMl  which  is  a hurizonral 
hve-asia  high-speed  milling  machine,  and  Figure  1-6  shows  the  UF-MTRC  HSM2 
vertical  three-osis  high-speed  milling  machine.  Both  of  (hese  machines  are  equipped  wilh 
open  archi lecture  control  systems  (hat  permit  Ihe  implcmemation  of  user-wriiien  servo 
algorithms.  The  HSM I serves  as  the  (es(  bed  for  this  study.  The  Y-axis  con.sists  of  a 
ballscrew  and  a rotary  motor.  This  drive  is  tuned  with  a FID  controller  with  digital  niters 
and  is  also  m-tuned  fora  sliding  mode  controller.  The  Z-axisof  the  machine  is  driven  by 
(WO  linear  motors  and  is  tuned  with  a simple  P)D.  Both  axes  are  used  lo  compare  and 
contrast  ballscrew  drives  and  linear  motor  drives. 


Figure  1-4.  Tracking  error 
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Figure  1-6.  High-speed  inilling  machine  (HSM2)  al  UF_MTRC 


CHAPTER  2 
LITERATURE  REVIEW 


Computer  Numerical  Cowiol  (CNC)  hui  become  ihe  most  imponam  pan  of 
modem  machine  tools.  With  the  continued  development  of  microprocessors.  CNC 
machines  have  become  more  technologically  advanced  in  their  computing  power.  Today’s 
machine  tool  controller!  may  possess  several  proces.sors:  some  for  generating  command.s. 
and  others  for  applying  control  algorithms  on  the  drives.  Some  controllcri  even  have 
dedicated  processors  for  every  drive  on  the  machine.  Machine  lool.s  have  entered  the 
digital  world,  and  with  this  new  technology,  advanced  control  stralegic.s  can  be  more 
easily  implemented  and  evaluated. 

In  order  to  test  control  stialcgies,  an  acceptable  cnntrollcrmusr  be  found. 
Conventional  CNC  technology  is  built  around  proprietary  urchitecture.s  with  incompatible 
software  and  operating  protocols.  Additionally,  these  controllers  are  expensive  and  costs 
loenhanceiheirflexibility  ore  prohibitively  high.  Open  architecture  controllers, 
generally,  allow  reseaichets  to  implement  und  evaluate  more  complex  control  strategics 
on  machines  in  less  time  and  at  a lower  cost  than  previously  possible.  With  Ihe  flexibility 
of  the  open  architecture  controller,  machine  tool  research  has  become  more  focused  on 
speed  and  precision. 

This  chapter  reviews  some  of  the  work  other  researchers  have  conducted  in  the 


control  of  machine  tool  drives. 


Conlral  Theory 


[n  generd,  two  types  of  comrot  strelegics  epplicd  TO  mechine  Tools-open  loop 
comrol  used  with  stepper  motors,  sod  closed  loop  (fcedbuck)  control  used  with 
servomotors.  Only  closed  loop  control  can  provide  the  degree  of  tracking  accuracy 
necessary  for  high*speed,  high'precision  metal  culling.  A closed  loop  feedback  controller 
compares  the  actual  output  from  the  drive  to  the  desired  input,  dclermines  the  error  and 
produces  a control  signal  to  eliminate  it.  [3|  The  most  widely  u.sed  controller  theory  i.s  the 
Proponionai.|megral*DerivaIive  IPID)  conuoiler. 

PID  Conirollcr 

The  basic  PID  controller  in  Figure  2*  I compares  the  po.sitional  feedback  to  the 
desired  position,  computes  the  error,  and  tries  to  eliminate  it.  This  controller  use.s 
proportional,  integral,  and  derivative  gains  to  adjust  the  performance  of  the  system. 

The  advantages  of  this  conlroiler  include  that  it  is  simple  to  install  and  it  is  ctusily 
understood  by  most  engineers. 

The  principal  disadvantage  of  a PID  controller  is  that  it  cannot  deal  with 
nonlinearilies.  Nonlitiearllies  which  arise  from  friction,  cutting  forces,  and  machine 
structural  modes  cannot  be  easily  compensated  by  a PID  controller.  These  nonlinearities 
may  cause  instability  which  llnuis  the  closed  loop  position  bandwidth  of  the  .system,  and 
results  in  low  performance. 

Fcedforwnrd  Ibchniques 

Some  researchers  have  modified  the  basic  PID  controller  to  compensate  for  some 
ofthc  nonlincaritcs  by  adding  feedforward  lerras  to  the  control  loop.  Tlusly  and  Chen  [4| 


dsall  wiUi  Oie  cancellaiion  of  vibralions  of  lowly  damped  slruclunil  modes  lhat  limii 
bardwidih-  They  suggesled  ihree  differem  ulicmaiives  lo  eliminaic  the  lowly  damped 
siniciural  modes.  The  firsl  alicmaiive  used  a feedforward  compeoMiion  filler  lhat 
cancelled  oui  ihe  lowly  damped  mode  and  was  iropicmenled  using  a roiaiy  encoder.  The 
second  aliemative  used  a llneai  encoder,  a similar  feedforward  compensation  filter,  and 
acceleromeiric  feedback  lo  increase  syslem  performance.  The  third  altemaiive  used  a 
tuned  passive  damper  lo  modify  the  dynamics  of  the  servodrivc.  The  di.sadvaniage  of  the 
first  and  ihird  ajiemalives  is  dial  ihey  ctmnoi  compen.sale  for  frequencies  of  llie  lowly 
damped  modes  ihal  change  wilh  ihc  raollon  of  the  machine.  In  other  words,  the 
feedforward  compensation  fillers  nre  tuned  to  cancel  only  one  specific  frequency. 

Week  15]  introduced  an  inverted  syslem  model  (IKF)  and  a low  pass  filter  to 
compensate  for  undesired  frequencies.  The  IKF  coniroller  can  increase  the  iracking 
performance  of  the  drivet  bui  it  is  only  useful  for  linear  sysiems.  Furthermore,  like  the 
TIusty  and  Oicn  model.il  cannot  readily  compensate  for  any  nonlinear  .struciurni  modes. 
Zero  Phase  Controller 

Tomizuka  |6)  inirodueed  the  zero  phase  error  tracking  controller  (ZPETC)  to 
improve  performance.  This  feedforward  coniroller  cancels  all  closed  loop  poles  and  all 
closed  loop  zeros  Ihal  are  inside  ihe  unit  circle.  It  also  cancels  all  phase  shifu  of  the 


eshibits  zero  phase  shift  between  the  desired  and  actual  output  at  all  frequencies.  Since 
this  controller  could  nolbcca.sily  implemented  on  an  nclual  servodrivc,  Tomizuka  and 
Yao  |7|  developed  an  adapiive  version  of  Ihe  ZPETC  Ihnl  uses  a pole  placement  controller 
nJgonlhm. 


Pole  Placement  Control 


Aliinios  and  Erkorkraoz  |8]  applied  a |>ole  placement  controller  shown  in 
Figure  2-2  with  Kalman  rillersond  friction  compensation.  The  values  of  the  polynomials 
WO.  J(0.  and  1\z)  were  determined  by  comparing  rhe  dosed  loop  transfer  function  in 
Equatton  2-1  to  a desired  model  B„{z)IA^z).  The  system  Sfof/lfal  was  identified  using 
white-noise  excitation  and  a least  squares  method. 

BT 

X,(e)  A^R*BS 

Kalman  filien  are  used  to  reject  noise  and  disturbances.  The  main  limitation  of  this  type 
of  controller  is  that  it  cannot  cope  with  any  changes  in  the  system  Blzi/Alz)-  For  example, 
if  the  system  is  not  precisely  Identified  or  is  nonlinear,  then  the  overall  system  may 
become  unstable. 

H-Infinlty  Control 

Van  Brussel  |9j  implemented  an  H-infinity  algorithm  to  control  linear  motors  with 
a varying  table  mass  as  a parametric  uncenainiy.  He  found  out  that  the  tracking  error 
becomes  too  high  when  the  mass  is  increased  to  300%.  He  also  concludes  that  the 
limitations  of  the  conuol  method  nre  due  to  poor  selection  of  weighting  functions. 

Sliding  Mode  Control 

Sliding  mode  control  theory  (SMC)  is  becoming  a favorite  for  mnehine  tool 
ie,searchers.  It  Is  a variable  slmciure  controller  that  does  not  require  prcci.se  system 
identification.  Sloline  |10]  describes  the  methodology  for  implementing  a sliding  mode 
controller,  and  rorcomecling  problems  like  chatter. 


Van  Brussel  [9]  implememed  an  SMC  and  found  ihai  ihe  tracking  error  remains  at 
6/rm.  even  if  the  moving  mass  on  a linear  moujr  was  increased  by  300%;  buldidnol 
discuss  the  overall  bandwidth  of  the  system.  Allimas  and  Erkorkmaz  |8|  used  the  SMC  to 
control  rotary  servomotors  and  were  able  to  achieve  a positional  bandwidth  of  100  Hz. 
Other  Control  Theories 

Boucher|ll|  used  a genemlizcdpredictivecitscade  control  (GPCC)  algomhm.  By 
predicting  Ihe  behaviorofthe  plant  for  a defined  bandwidth,  he  generated  a set  of  future 
control  values  to  minimize  thcemor.  Keren  and  Lo  [12]  used  a variable.gain  cross 
coupling  controller  that  generated  correction  signals  for  each  axis  based  on  contounng 
emots  supplied  through  feedback  from  all  axes  on  the  machine. 

Linear  Motors 

Some  of  the  conventional  drives  on  high-speed  milling  machines  are  being 
substituted  by  direct  drives  that  use  linear  motors.  Prilschow  |I3]  discussed  the 
advantages  of  linear  motors  over  conventional  drives  in  high-speed  machining,  as  well  as 
the  dincrences  among  the  different  types  of  linear  motors.  He  concluded  that  because  of 
Ihe  absence  of  compliant  transmission  elements,  linear  motors  can  have  larger  position 
loop  bandwidiha  and  higher  accelerations  than  conventional  rotary  servomotors. 
Furdtermore,  linear  motors  can  deliver  higher  speeds  which  increase  the  productivity  of 
Ihe  machine.  Finally,  he  notes  that  DC  linear  molars  are  more  suitable  in  machine  tool 
□pplication.sihan  AC  induction  type  linear  motors  because  of  their  ability  to  deliver  a 
higher  force  per  unit  weight  of  Ihe  magnets  and  coils. 


CHAPTERS 

DESIGN  FOR  LINEAR  MOTION 

Linear  motion  in  high  perrormoncc  milling  mochines  cenioni  on  using  two  basic 
axis  drive  designs:  the  leadscrew  drive  and  Ihe  linear  moior  drive.  Thcleadscrew  drive 
olTers  a low-cost  drive  with  high  slifTnikss;  and  because  ollhe  transmission  ratio,  it  is 
insensitive  to  changes  in  woritpiece  mass.  The  linear  motor  drive  offers  higher  spcecU  and 
higher  accelerations,  but  since  it  is  a direct  drive,  it  is  very  sensitive  to  changes  in 
workpiece  mass  and  to  external  disturbances.  In  this  chapter,  both  drives  are  discussed  to 
evaluate  their  strengths  and  weaknesses. 

Leadscrew  Drives 

A leadscrew  drive  is  the  most  commonly  used  actuator  for  high  performance  axis 
drives  in  milling  machines.  The  drive  in  Figure  T-l  shows  a leadscrew  drive.  The 
leadscrew  is  coupled  to  a rotary  .servomotor  that  rotates  the  leadscrew  in  Ihe  bearings.  The 
(able  and  nut  are  driven  down  the  helical  thread  of  the  leodscmw. 

Components  of  Leadscrew  Drives 

Direcl-curreni  rotary  servomotors  ore  the  most  widely  used  type  of  motors  m 
milling  machines.  Direcl-curreni  servomoiora  arc  better  suited  for  machine  tool 
applicaUons  than  AC  motors  because  they  have  a lower  moment  of  inertia  and  can  provide 
higher  peak  accelerations  than  AC  servomotors. 
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The  servoinoior  is  eonneeled  lo  Ihe  leadscrew  using  u flexible  disk  coupling  ihui 
allows  for  small  angular  and  axial  misalignmenl  between  ihe  molor  shaft  and  leadscrew 
end  yel  provides  high  lorslonal  rigidily.  Unfonunalely,  to  obuiin  the  high  torsional 
stifTness,  a large  coupling  is  used  which  has  a high  moment  of  inema 

The  leadscrew  is  held  in  two  seta  of  end  bearings.  The  bearings  - typically,  a pair 
of  angular  contact  ball  bearings  or  tapered  roller  bearings -closest  totlie  motor  handle  all 
of  the  axial  loads  imposed  on  the  leadscrew.  The  "free"  end  of  the  leadscrew  is  held  by  a 
radial  bearing.  The  radial  bearing  allows  for  thermal  expansion  of  the  leadscrew  and 
raises  the  minimum  critical  speed  of  the  leadscrew. 

Conventionally,  the  ballscrew  is  the  preferred  leadscrew  in  high-speed  drives. 
Precision  ballscrews  have  three  components;  a screw,  a nut.  and  recirculating  balls.  The 
balls  recirculate  throagh  concave  helical  grooves  ground  in  both  the  screw  and  the  nut. 
Backlash  is  eliminated  by  squeezing  Ihe  balls  between  Ihe  screw  and  the  nut.  This  is 
achieved  by  having  a nm  with  two  halves  which  have  slightly  mismatched  grooves.  The 
ballscrew  is  limited  in  speed  and  acceleration  because  of  the  recirculming  ball.s. 

An  alternative  to  the  ballscrew  is  the  hydroslatic  leadscrew.  High-pressure  oil  is 
supplied  10  pockets  that  are  machined  in  the  threads  of  the  nut.  This  oil  suspends  the  nut 
from  the  screw  which  make.sihe  hydrostatic  drive  non-tomact  and  nearly  friclionicss. 
The  speed  of  the  hydrostatic  leadscrew  is  limited  by  hydrostatic  fluid  cavitation  or  pocket 


Advantages  and  Disadvantages  of  Leadscrew  Drives 

The  Ann  disadvantage  is  the  limited  life  of  a ballscrew  due  lo  wear  caused  by 
ntechanical  contact.  This  is  not  a problem  for  the  hydrostatic  leadscrew. 
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The  second  disadvantage  is  limited  speed  and  travel  based  an  the  screw  dtameier 
and  length.  A ballscrew  is  limited  by  critical  speed  or  ball  speed.  A longer  Icadserew  of 
fised  diameter  has  a lower  critical  speed  depending  on  end  bearing  arrangetttent.  A larger 
diameter  screw  has  a lower  max  imam  operating  speed  because  of  the  higher  peripheral 
velocity  of  the  balls. 

A third  disadvantage  is  the  limited  position  bandwidth  that  can  be  achieved  with  a 
leadscrew.  Leadscrews  can  have  tontiona]  vibrations  which  limit  the  bandwidth.  Low 
pass  niters  are  commonly  used  in  machine  tool  drives  to  attenuate  higher  frequencies  and 
stabilize  the  control  loop.  However,  low  pass  niters  were  not  used  by  Tohme  |I4|  and 
Smith  |IS|  because  they  decrease  the  performance  of  the  drive  by  decreasing  the  overall 
bandwidili  of  the  system.  The  peaks  were  eliminated  by  in.seningan  inverse  filler  in  the 
feedback  loop  that  canceled  both  resonances.  These  digital  filters  extended  the  position 
loop  bandwidth  up  to  100  Hz.  The  problem  with  such  a design  is  that  a change  in  the 
drive  dynamics  can  allow  a mismatch  between  the  filierv  and  the  open  loop  response  Ihu.s 
causing  an  unstable  system.  Also  an  additional  controller  was  required  to  handle  the 
calculations  of  the  niters  for  every  servo  loop. 

A fourth  disadvantage  is  the  location  of  the  feedback  devices.  From  Figure  2-1 
rotary  encoders  are  typically  mounted  on  the  backend  of  rotary  motors  to  provide  position 
feedback.  The  rotary  encoders  may  not  sense  the  actual  position  of  the  table  due  to 
mechanical  vibrations  outside  the  control  loop  or  thermal  expansion.  This  can  lead  to 
inaccuracies  in  the  machined  port.  An  allemalive  is  mounting  linear  encoders  parallel  to 
the  motion  of  the  table  that  can  report  accurate  table  positions.  The  linear  encoder  senses 
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vibrations  which  were  pieviously  outside  the  position  loop.  These  vibrations  may 
destobilizc  the  position  loop  if  uncompensated. 

A disadvantage  of  a hydrostatic  screw  is  that  it  requires  u continuous  supply  of 
cool,  high  pressure  oil.  This  requires  u supply  pump,  an  oil  reservoir,  a heat  exchanger 
and  3 chiller  Additionally,  the  hydrostatic  oil  must  be  captured  upon  leaving  the  nut  and 
must  be  isolated  from  other  fluids. 

Advantages  of  a Icadsctcw  drive  includes  low  cost,  simple  assembly  procedures, 
low  sensitivity  to  changc.s  in  workpiece  mass,  and  high  static  and  dynamic  siufTiness. 

Linear  Motor  Drives 

A linear  motor  drive  Is  a direct  drive  which  acts  directly  on  the  tabic  mass  without 
a transmission.  Linear  motors  are  mounted  parallel  to  the  guideway.s  with  the  shorter 
motor  assembly  attached  to  the  table  and  the  longer  assembly  attached  to  the  bed.  An 
optical  linear  encoder  is  mounted  parallel  to  the  motion  and  provides  position  feedback. 
The  direct  drive  eliminates  oil  mechanical  transmissions;  therefore,  all  the  performance  of 
the  drive  is  limited  to  one  component  - the  linear  motor. 

Fundamentally,  linear  motors  ate  rotary  motors  rolled  out  flat  as  shown  in 
Figure  3-2.  The  moving  part  is  typically  called  the  rotor  and  the  fixed  part  the  stator. 
Linear  motors  ore  divided  into  two  main  groups;  AC  motors  and  DC  motors.  AC  motors 
are  also  divided  into  two  subcatngories:  the  linear  inductance  motor  and  the  linear 
reluctance  motor  Generally,  AC  motors  are  used  for  movement  over  long  distance.s  as  in 
rail  locomotives  while  DC  motors  are  used  for  precise  movement  overshonet  distances. 
DC  moiois  ore  divided  into  three  subcategories:  the  moving  magnet,  ihe  iron  core 
irontess  motors. 


e,  and  the 


CumponcnU  of  DC  l.intar  Motor  Drivci 

The  iwo  pnoiary  coraponenis  of  brushless  DC  linear  moiois  are  ihe  permanem 
magiiei  assembly  and  the  coil  assembly.  The  basic  theory  behind  the  operation  of  linear 
moioR  is  similarto  that  of  rotary  motors.  When  the  conductors  of  n ccnain  leitgih  f in  the 
coil  carries  a rrurreni  / in  a direction  perpendicular  to  the  iniigneiic  field  of  flu*  density  B. 
ihecoil  experiences  a force  Fpeipendicular  to  both  the  direction  of  the  current  flow  and 
the  magnetic  Held.  Figure  3-3  shows  the  orientation  of/,  B.  andF.  The  magnetic  field  Is 
generated  by  permanent  magnets.  The  magnets  used  on  the  linear  motors  at  UF-MTRC 
high-speed  milling  machine  are  a neodymium-iron-boron  ceramic  magnets  which 
produces  a stronger  field  than  sumdatd  magnets. 

From  Basak  (161,  the  bnishlcis  DC  linear  motors  are  divided  into  three  categories; 
moving  magnet  motors,  moving  armature  motors  (or  iron  core  motors],  and  moving  coil 
motors  (or  ironless  motors). 

A moving  magnet  motor  ia  similar  to  the  bru^less  rotary  motor.  Figure  3-4a 
shows  how  the  stator  is  made  from  laminations  while  ihe  raiof  is  made  from  permanent 
magnets.  These  motors  have  a high  powerless  making  them  very  inefficient  compared  to 
the  moving  armature  motor. 

The  moving  armature  motor  of  Figure  3-4b  has  stationary  magnets  and  moving 
Inminaiions.  They  are  more  efficient  than  the  moving  magnet  motors  because  the 
windings  are  opposed  by  magnets  at  nil  lime. 


The  third  type  is  the  moving  coil  motor  of  Figure  3-dc.  It  is  similar  to  the  moving 
armature  magnei  escepi  Ihe  laminaiions  are  free  of  iron  which  eliminates  cogging  force. 
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Magnets  on  both  sides  of  the  motor  ore  used  to  generate  more  thrtisl  foroe  and  provide  a 
flux  return  path  for  the  magnetic  cireuit. 

For  this  work  the  motor  that  best  fits  the  application  and  constraints  is  the  iron  cure 
motor.  Tobic  3*1  shows  some  of  the  differences  between  the  three  motors.  Aiihoughthe 
charocicrtstics  of  the  iron  less  motor  are  desirable,  the  thrust  force  of  the  available  siaes  is 


loo  weak  for  the  application.  Because  of  the  inefficiency  of  the  moving  magnet  motor,  the 


iron  com  motor  is  the  best  suited  motor. 

Table  3-1.  Comparison  of  brushless  direcl-currenl  linear  motors. 


Moving  magnet  Iron  core 
Weight  of  moving  unit  Medium  High 

Cost  Moderate  Low 

Cogging  force  Prc.vent  Present 

Thnislforcecapabiliiies  High  High 

EfTiciency  Low  Medium 


Iron less 
High 


High 


Advantages  and  Disadvantages  of  Linear  .Motor  Drives 

One  disadvantage  of  large  linear  motors  is  heal  generation.  The  coils  ore  veiy 
large  in  size  and  need  to  be  cooled  in  order  to  achieve  the  highest  efricicney.  The  required 
cooling  system  increases  the  cost  and  complexity  of  the  machine.  Also  temperature 
sensors  are  needed  to  monitor  the  motor  continuously. 

The  second  disadvantage  is  that  linear  motors  have  low  efficiencies.  This  low 
efficiency  is  mainly  due  to  large  air  gaps  between  the  sialor  and  the  loior.  This  gap  allows 
a high  flux  leakage.  From  Basok  [16],  the  emdency  of  some  linear  motor  could  be  lcs.s 


than  30%. 


A third  disodvuruage  U the  high  lateral  attractive  fotte.  The  i 


unlike  [he  ironle&s  motors  have  a large  attraction  between  the  magnet  and  the  steel 
laminations  in  the  coil.  These  foircscan  be  as  much  o.s  three  limes  largenhan  iheihruM 
force.  The  attractive  force  of  the  moiom  selected  for  this  work  i.s  aboul  35.000  N while  the 
thrust  force  is  only  I3.000N.  Thisaitmciive  force  can  esert  a large  force  onto  the 
guideways  increasing  friction  and  wear  in  the  guideways.  It  is  possihlc  in  cancel  this 
allraciive  force  by  mourning  two  linear  motors  back  lo  back.  This  large  force  also  causes 
problems  during  installation.  Mounting  procedures  ntusl  prevent  the  coils  and  the 
magnets  from  touching.  If  the  coil  gels  stuck  on  the  magnets,  removal  may  be  impossible 
without  damaging  the  motor. 

The  fourth  disadvantage  is  that  linear  motors  must  be  protected  from  chips  ihai  are 
ulimcted  to  the  mogrreiic  components.  Metal  chips  can  be  pulled  into  the  air  gap  leading 
to  a short  and  damaging  the  coils. 

The  fifth  disadvantage  in  linear  motors  is  that  the  lublc  cannot  be  held  in  position 
without  power.  Them  is  nothing  loslop  the  motion  of  the  mass  in  case  of  a loss  in  power 
Additionally,  an  cMemal  brake  is  required  for  vertical  uppheattons. 

The  sixth  disadvantage  is  cogging  forces  in  the  iron  core  motor.  These  cogging 
forces  are  ahout  5%  of  ils  conlinuous  force.  These  disturbances  cause  unaccepitUjIe 
position  error  and  must  be  compensated. 

The  seventh  disadvaniage  is  the  high  sensitivity  to  changes  In  workpiece  muss  or 
exiemol  force  disrurbonces  like  culling  forces  and  friciion  forces.  Figure  3*5  shows  the 
step  response  of  a linear  moior  for  two  di^ercni  musses  without  retuning  the  servo  loop. 
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ThetmusinaiscSilS  twice  chemise  in  case  A.  Siniiiarly.  external  forces  will  cause 
position  changes. 


The  main  advantage  of  diieci  drives  Is  the  speed  and  acceleration  they  can  provide. 
Linear  motors  can  achieve  .speeds  up  to  2S0  m/min,  about  5 times  mom  than  IcotLscrcw 
drives.  They  can  achieve  high  accelerations:  up  to  10  g unloaded.  High  bandwidth  is 
easily  achievable  with  ininitnal  control  work. 

Another  major  advantage  is  the  slimness.  Overall  sti^ness  is  highly  desirable  in 
machine  tools.  Reportedly,  linear  motors  can  have  higher  stiffnesses  than  lead.screw 
drives  [13].  Stiffness  in  linear  motois  is  dependent  on  the  position  loop  tuning.  High 
gains  In  the  position  loop  can  increase  stiffness  in  the  drive. 

Sunimory 

A comparison  between  a linear  drive  and  a leadscrew  drive  is  shown  in  Table  3-2. 
In  summary  the  linear  motor  drive  can  achieve  a higher  acceleration  than  a ballscrew 
drive,  its  major  drawback  is  its  sensitivity  to  mass  changes  and  external  disturbances. 


Table  3>2.  Comparison  between  linear  molor  drive  and  leadscrew  drive 


Linear  motor  drive 


High  speeds  up  [o  2SD  m/min 


Acceleration 

Travel 

Sensitivity  to  mass 
changes 

Sensitivity  to  ester- 
Cogging  force 

Backlash 
Dynamic  stIITness 
Installation 
Operating  life 
Efncieney 


High  (Limited  by  the  motor 
peak  occclcraiian) 
Unlimited  travel 
High 

High 

High  (t^ut  5%  continuous 

High  (Depending  on  tuning) 

Complicated 

High  (No  wear) 

High 


Leadscrew  drive  using  a brush- 
less rotary  motor 
Limited  to  the  leadscrew's 
maximum  speed 
Limited  by  flexibilities  of  the 

Limited  travels 


Moderate  (about  l%continu- 

High 
Simple 
Limited  life 
High 


CHAFFER  4 

TUNING  ROTARY  MOTOR  DRIVES  FOR  WIDE  BANDWIDTH 


The  geomeiry  ofmosc  aircmR  pons  includes  deep  pockets  with  small  dimensions. 
High-speed  machining  these  pockets  lequiies  high  accelerations  in  ordenhai  the 
maximum  speed  con  be  reached  over  Ihe  short  distance.  A common  aircraft  pan  requires 
more  than  SOO  acceleration  and  deceleration  commands  in  less  than  one  minute.  Only  a 
feed  drive  with  wide  enough  bandwidth  can  accurately  track  such  commands.  Wide 
bandwidth  is  not  easily  achievable  because  structuiaJ  modes  in  the  dnvc  limit  bandwidth. 

As  rrrenlioned  earlier,  there  are  two  methods  to  achieve  a high  bandwidth  feed 
drive:  altenng  or  improving  the  control  method  to  deal  with  these  structural  modes,  and 
redesigning  the  feed  drive  axis  using  linear  motors.  The  next  chapter  deals  with  Ihe 
coninal  of  a linear  motor,  while  this  chapter  investigates  the  source  of  the  structural  modes 
and  methods  to  treat  them  In  the  control  system. 


The  high-speed  milling  machine  HSM I at  the  MTRC  is  used  as  a test  bed  to  show 
the  differences  between  the  two  different  drives.  The  machine's  Y-axis  illustrates  the  use 
of  alradllionni  ballscrew  drive.  The  feed  drive  consists  of  both  elecirical  and  mechanical 


systems  which  include:  the  controller,  PWM  amplifier  and  power  supply.  Ihe  servomotor, 
encoder,  tachometer,  ballscrew,  ballnul  bearings  and  the  moving  mass. 


AnnJysis  of  l^adwrew  Kccd  Drives 


The  ballscrew  drive  is  shown  in  Figure 4- 1.  The  ballscrew  is  driven  by  a 
pcrmimcnl  magnet  brushless  DC  servomotor  made  by  Kollraorgen,  model  number  B- 
804B.  The  servomotor  has  permanent  neodyraium-iron-boion  alloy  magnets  ttiiached  to 
theiotary  shaft  and  a six-pole,  three-phase.  Y-stator  windings  attached  inside  a cast 
aluminum  housing  A resolver,  a tachometer  and  an  encoder  are  mounted  at  one  end  of 
the  motor  shaft,  and  on  the  opposite  side  the  ballscrew  is  coupled  through  a flexible 
coupling.  The  motor  is  driven  by  a PWM  (Pulse  Width  Modulation)  ampliricr  which  is  set 
to  operate  in  the  current  mode.  The  resolver  is  fed  bach  into  the  amplifier  and  is  used  to 
provide  a commulalion  signal.  The  velocity  loop  is  closed  with  a 17  V/krpm  tuchometer 
and  the  position  loop  is  closed  through  an  A-B  quadrature  encoder  with  20,000  counts  per 
revolution.  Both  the  velocity  and  position  loops  arcclosed  in  the  controller.  Theflexible 
coupling  is  altached  to  a ballscrew  that  is  S2  inches  long.  The  ballscrew  is  ailached  to  the 
driven  load  or  table  through  a baJlnut.  The  lead  of  the  ballscrew  is  20  mm  per  revolution. 

The  controller  uses  from  two  separate  components.  The  first  component  is  the 
commercial  CNC  conlroller  PM  AC  by  Delia  Tau.  This  controller  is  used  for  data 
acquisition,  command  generation,  and  other  CNC  housekeeping  function.s  like  executing 
PLC's  and  event  moniloring.  This  cottlroller  uses  the  Motorola  S6300  fixed  point  DSP 
chip,  which  is  not  capable  of  handling  the  time  processing  of  high  order  digital  fillers. 

The  .second  controller  is  a commercial  floating  point  DSP  board  which  is  used  only  for  the 
control  loops.  This  DSP  board  allows  the  implementation  of  any  control  strategy  as  long 
as  its  compulation  time  for  nil  ihe  axes  does  not  exceed  the  servoloop  closing  lime.  The 
servoloop  lime  is  sot  at  22 1 ps. 
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Servomotor 

SD 


Figure4-2.  Block  diagrjini  of  the  opcn*lGop 


AMfLlHER 


Figure  4<3.  Block  Diagram  of  irarirrrerfuncLjotr 


Table  4-1  shows  The  values  of  conslams  in  Figure  4-3  und  Equaiion  4-2  os  supplied 


by  The  monufaciurcrs. 

Table  4-1,  Values  for  ballscrew  drive  porameicrs 


Ajopliner  uuioff  frequency 
Tachometer  low  pass  flhcr  frequency 
Ampliner  gain 

Tachometer  sensillviiy  gain 


®co  ~ ^300  rad/s 
«„=  !2600rad/s 
*■0=  15.1  AA< 
Kr=  1.546  Nm/A 


//7=  0.0477  V/(nid/s) 
J = 0.00840  kg.m- 


The  Transfer  fuiTclion  of  (he  actual  syslem  is  then  obuined  by  performing  a sine 
sweep  (esT  on  The  feed  drive.  A sine  sweep  lest  is  a method  used  to  determine  The  transfer 
function  of  a system.  A series  of  sine  waves  with  equal  amplitudes  and  varying 
frequencies  are  generated,  stored  on  the  controller,  itnd  then  setil  to  the  servomotor  via  the 
amplifier  in  an  ascending  order.  As  the  wave  causes  the  motor  to  oscillate,  the  feedback 
sensors  ore  sampled  by  the  controller.  In  ihisca.se.  the  feedback  sensors  are  the 
tachometer  ond  the  encoder.  After  all  the  dttla  is  collected,  u Fourier  transform  in  applied 
to  the  sampled  tachometer  response  und  the  result  of  the  output  over  the  input  is  the 
transfer  function. 

The  result  of  the  sine  sweep  test  is  the  transfer  function  plot  in  Figure  4-4. 

The  figures  of  the  transfer  functions  show  twoseparaie  plots:  the  magnitude  and  the 
pha.se.  The  bandwidth  ofasystem  is  defined  us  the  frequency  range  before  the  magnitude 
drops  below  -3db.  The  phase  plot  shows  the  phase  shift  at  each  frequency.  The  system  in 
Figure  4-4  has  u slope  bccau.se  this  is  Ihc  result  us  seen  from  the  tachometer,  which 
measures  velocity.  The  encoder  is  also  sampled,  but  because  die  tachameier  Is  more 


sensiUve,  il  is  used  to  generate  die  tiansrer  Functions.  Integration  of  the  loehomctcr  sigruti 
yields  die  position  signal,  llie  figure  .shows  that  there  arc  no  noticcabie  vibrations  under 
i.OOO  Ha.  ihecutofT  fiupuency  of  the  ampilfier.  This  means  that  up  to  the  cutoff  Frequency 
the  relationship  between  the  command  and  the  oetuai  vciocity  is  one  to  one.  A 3,000  rpm 
command  yieids  a 3,000  rpm  response.  At  Frequencies  higher  than  the  cutoff  frequency 
the  mlutionship  between  the  command  and  acliiai  veiocity  is  attenuated.  The  phase  piot 
represents  the  time  shift  ora  time  deiay  between  the  command  and  actual  velocity. 

After  the  sine  sweep  lest  is  completed,  the  ballscrcw  is  coupled  to  Ihe  motor.  A 
second  sine  sweep  test  is  esecuted  and  Figuiv  4-5  shows  the  open  loop  transfer  function  of 
the  system  with  the  ballscrcw.  As  expected  several  dynamics  associated  with  the  coupled 
load  appeared.  These  plots  show  modes  at  50,  550  and  1,050  Hz.  These  modes  cause 
insiobiiiiy  in  the  control  loop  If  a 100  Hzbaadwidth  is  to  be  achieved.  Removal  of  them 
requires  identifying  their. sources. 

Lcadscrcw  Modal  IdenUncation 

In  order  to  understand  the  behavior  of  the  system,  a more  thorough  understanding 
of  the  system  is  required  which  leads  to  the  developmem  of  a more  accurate  transfer 
function  model.  The  modes  in  Figure  4-5  arc  from  lorsional  and  axial  flexibilities 
associated  with  Ihe  table,  ballnui.  ballscrcw,  coupling,  and  motor  assembly.  Tliese  lowly 
damped,  second-order  vibrations  can  cause  the  system  to  become  unstable  when  they  are 
introduced  in  to  the  closed  position  loop.  In  most  conventional  machine  tool  applications, 
low-pass  filters  are  inserted  into  Ihe  position  loop  to  auenuaie  these  vibrations,  hut  they 
also  limit  the  system's  bandwidth.  The  objective  in  this  work  is  to  design  digual  fillers 
that  remove  Ihesc  vibrations  and  still  keep  a wide  bandwidth  up  to  100  Hz.  First  a model 
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of  these  second  order  dynmnics  imisi  be  added  to  the  discrete  ininsler  ruitcilon  in 
Equation  4-2.  Figure  4-6  shows  the  free  body  diogram  of  a second  order  system  model. 
From  Newton’s  second  law  of  motion,  Equntion  4-3  is  obtained 


mX*c(X-*-)  + iK;f-;if,)=  0 
A Laplace  ininsform  to  Equation  4-3  yields  the  following. 

Xfr)  cs*k 
^r(»)  ' a*  k 


(4-3) 


(4-4) 


m = The  mass  of  the  column  (kg), 
r = Damping  coefRcienl  (N/(m/sl). 
k = Spring  stiffness  (N/m). 

Each  of  the  three  modes  is  modeled  as  a second-order  system  as  shown  in 
Equation  4-4.  They  are  then  cascaded  together  to  form  the  entire  transfer  function.  These 
vibrations  ore  then  classified  into  two  types:  inside  the  loop  and  outside  the  loop. 
Vibrations  inside  the  loop  cause  system  instability;  these  have  to  be  removed  in  order  to 
achieve  wide  bandwidth.  Vibrations  outside  the  loop  do  not  cause  .stability  problems,  but 
they  might  leave  undesirable  machining  marks  on  the  ports.  Figure  4-7  shows  the  system 
G^s)  with  two  second  order  modes  that  are  added  to  the  inertia  of  the  system  and  are 
considered  Inside  the  loop.  The  vonables  in  Figure  4-7  are  derived  from  the  variables 
from  Equation  4-4.  Equation  4-5  shows  the  relationship.  The  values  of  these  variables 
are  computed  from  measurements  conducted  on  the  actual  system  and  from  data  supplied 
by  manufacturers.  The  value  of  Ihc  total  system  moment  of  inertia  is  caiculaied  by 
summing  all  of  the  inertia  of  the  moving  parts  in  the  system. 
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= Natural  frequency  (rad/s). 

^ = Damping  ratio  (N/m). 

This  includes  the  moments  of  Inertia  of  the  following:  the  rotor  or  motor  shaft,  coupling, 
ballscrew,  ballnut.  and  table.  The  table  and  ballnut  inertia  onr  calculated  as  reflected  on 
the  motor  shaft.  Table  4-2  shows  the  moments  of  inertia  of  alt  the  components  of  the  Y- 
axis  on  HSMI. 

Table  4-2.  Moments  of  inertia  of  the  ballscrew  drive 
Description  Moment  of  inertia  J (Kg.m^)  Percentage 

BnlUciew  0.00182  S6.6 

Rolor  0.00340  26 

Table  and  nut  0.00820  2.5 

Coupling  0.00485  15 

Total  0.03227  100 

Figure  4-8  shows  the  open  loop  model  with  C/s)  that  includes  two  structural  modes 
added  to  it.  The  remaining  values  are  then  extracted  from  actual  measurements. 

Two  of  the  modes  that  are  in.side  the  loop  are  suspected  to  be  torsional  flexibilities 
in  the  ballscrew  and  motor  assembly.  A test  is  pertbrmed  to  properly  identify  them  Two 
accelerometers  mounted  opposite  each  other  are  used  to  sense  torsional  vibrations  and 
reject  axini  vibrations.  The  same  sine  sweep  test  used  for  the  previous  open  loop 
measurements  is  conducied  again.  The  lesult  is  shown  in  Figure  4*9  which  shows  Ihc 


mode  shapes  of  the  ballscrew  at  Ihe  550  und  1 , 1 50  Hz  modes.  This  test  proves  i 
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Figure4-6.  Spring'MaM-Damper model. 


aiii»  1 bL» 

Figure  4-7.  Dynumicjf 


AMPPFIEK 


Figure  4<8.  Open  loop  with  dynamics 


(dashed  line). 


niKr  is  ctiosen  for  removing  ihc  undesirable  siniclurol  modes  from  ihe  ironsfer  function. 
A short  description  of  Ihe  fillers  is  given  in  Appendin  B.. 

Table  4-3.  Parameters  evaluated  from  experimenul  measurements 


k,  3.80e5  Nm/rad 

CD„i  3220rad/s 

0.02 

k2  I.52e6  Nm/rad 

(0;,2  6370  rad/s 


Stiffness  of  mode  I 
Natural  frequency  of  mode  1 
Damping  ratio  of  mode  1 
Sliffnctts  of  mode  2 
Natural  frequency  of  mode  2 
Damping  ratio  of  mode  2 


Digital  fillers  arc  needed  to  atienuule  the  siruciurai  modes.  Digitui  filters  can  be 
miroduecd  anywhere  in  the  digital  part  of  the  control  loop.  Three  differenl  filters  were 
chosen  to  Improve  the  feed  drive  performance.  Figure  4-1 1,  shows  two  of  these  positions, 
the  velocity  loop  filler  and  the  position  loop  filter.  The  velocity  loop  filler  W^(z)  is 
installed  inside  the  velocity  loop  and  is  used  to  filter  Ihe  Iwo  main  torsional  modes.  The 
position  loop  filler  ff„U)  is  installed  outside  the  velocity  loop  and  inside  the  position  loop, 
its  purpose  is  to  clean  anything  lefl  untreated  in  the  closed  velocity  loop.  A third  tiller  is 
placed  outside  Ihe  position  loop.  That  type  of  filtering  is  called  feedforward  nilering  since 
it  does  not  affect  the  feedback  loop.  Feedforward  can  be  used  to  improve  on  the  system  by 
altering  the  commanded  motion  and  it  does  not  affect  stability.  The  implementation  of 
feedforward  is  described  in  Chapter?. 

The  objective  of  the  digital  fillers  is  to  cancel  the  two  torsional  nexibililies.  Each 
mode  is  represented  by  Iwo  poles  and  two  zeros  near  the  unit  circle  in  the  z-plane.  The 


digital  filter  ihccreiically  uses  poles  lo  cancel  structural  mode  zeros.  Likewise,  the  filler 


uses  zeros  to  cojiccl  strtictuRl  mode  poles.  The  mult  should  equal  one.  ll  is  possible  to 
design  an  unstable  filler  If  die  Tiller  poles  are  matched  to  zeros  ihai  are  outside  the  unit 
circle.  Consequenil)'.  the  poles  am  placed  as  close  as  possible  to  the  zeros  but  in  the  unit 
circle.  The  repmeniation  of  the  digital  Tiller  is  described  in  the  following  equation. 


I -|2r.ieose,,|z-'  + r.|V’ 
l-|2r^,cos0,„|z-l+r„V^ 


(4-7) 


The  poles  and  zeros  are  described  by  r the  radius  from  the  center  and  0 ihe  angle, 
is  ihe  transfer  function  for  one  sets  of  poles  and  zeros. 

Table  4-4  shows  how  the  poles  of  Ihe  syslem  are  the  zeros  of  the  Tiller  and  how  the 
zeros  are  now  Ihe  poles.  The  same  open-loop  sine  sweep  test  is  performed  with  Ihe  niters 
implcmenied.  Figure  4-12  show.s  the  compensated  open-loop  transfer  function.  Both 
peaks  from  the  pievious  measurement  are  now  rcmovedi  therefore  Ihe  velocity  loop  can 
be  closed  wiih  Ihe  highest  possible  gam.  Once  the  velocity  loop  is  tuned  to  the  best 
achievable  gain  A',  a uansfer  function  Is  measured  by  applying  a. similar  sine  sweep  test. 
Figure  4-13  shows  the  compensated  velocity  loop  of  the  system.  The  bandwidth  of  the 
velocity  loop  Is  around  2IX)  Hz.  with  a phase  shift  of  l(K)  degrees  at  200  Hz.  The  50  Hz 
mode  remains  although  it  does  not  affect  the  stability  of  the  system.  Several  attempts  to 
remove  it  in  Ihe  position  loop  filter  7/p(c)  failed  because  the  mode  is  non-linear.  As  the 
lable  traverses  the  rails,  the  natural  frequency  of  the  that  mode  changes. 

The  position  loop  is  closed  without  any  Tillers  added  to  it.  Figure  4-14  shows  the 
transfer  function  of  the  position  loop  between  Ihe  input  and  the  encoder  as  an  output.  The 
position  loop  reached  around  100  Hz  although  the  notch  due  to  the  50  Hz  mode  did 
deteriorate  Ihe  performance.  The  1 00  Hz  is  achieved  with  Ihe  highest  possible  position 
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loop  gain.  The  gains  are  reduced  in  ihe  final  impiementaiion  as  a .safeguard  against 
unforeseen  changes  in  the  machine  performance  producing  instabiiity.  The  current 


bandwidth  after  this  safeguard  is  70  Hz. 
Table  4-4.  Compensation  poles  and  zeros 


Poles  and  zeros  of  C^l'. 


Poles  and  zeros  of  W,,(z) 


Radius  Angle  Radius  Angle 


Angle  Radius  Angle 


0.98588  +0.7108  0.97319  +0.5101 

0.98588  -0.7108  0,97319  -0,5101 

0.97139  +1-4508  0.95871  +1.2329 

0.97139  -1.4508  0.95871  -U329 


0.97319  *0.5101  0.98588  +0.7108 

0.97319  -0.5101  0.98588  -0.7108 

0.95871  +1.2329  0.97139  +1.4508 

0-95871  -1.2329  0.97139  -1.4508 


After  the  position  loop  Is  closed  some  motion  commands  are  performed  to  show 
Ihe  performance  of  Ihe  machine  Figure  4-15  shows  a 2 inch  (50800  cts)  move  at  a 
feedrate  of  200  inymin  and  a 1 g acceleration.  Figure  4-16  shows  the  sieudy  state  error  of 
the  system.  The  steady  slate  error  could  be  better  seen  in  the  following  plot.  Figure  4-17. 
which  shews  that  Ihe  following  error  is  about  500  counts  (0.02  inches).  This  error  only 
shows  dial  die  actual  position  is  lagging  behind  the  commanded  position  by  0.02  inches. 

Ti  is  equivaleni  to  a lime  delay.  Since  il  is  a constant  during  steady  state  motion  it  should 
not  affect  the  accuracy  of  the  pan.  This  lag  error  is  only  significant  on  machined  puns 
when  several  axes  are  commanded  to  make  u coordinated  move.  For  instance,  if  each  axes 
has  a different  lime  delay  between  the  command  and  Ihe  actual  position,  a circular  move 
will  be  actually  elliptical.  Tuning  all  axes  to  have  similar  following  eirors  is  not  possible 
since  they  have  different  components.  Feedforward  is  later  applied  to  reduce  these  errors. 
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Figure  4-1 5.  Commanded  and  eciuol  morion.r 
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I(  is  evidenl  lhal  with  careful  tuning  and  good  undersuuiding  of  a system  a wide 
bandwidth  is  possible  on  ballscrew  drives.  The  use  of  diglutl  niters  has  iticreiised  the 
bandwidth  of  a ballscrew  drive  from  30  Hz  !o  100  Hz.  This  more  than  triples  the 


CHAPTER  5 

HIGH  BANDWIDTH  USING  SLIDING  MODE  CONTROL 

Tic  digital  niters  insened  in  ihe  PID  loop  suffer  from  ihc  lack  of  robustness.  The 
niiers  were  designed  loesaclly  Til  die  measured  irunsfer  funcUon.  If  ihereare  any  changes 
in  Ihe  transfer  funclion  the  fillers  will  slop  cancelling  Ihe  poles  and  aeroes  and  ihcy  will 
cause  an  insiablliiy.  A sliding  mode  coniroller  is  a robust  controller  that  should  provide 
high  bandwidth. 

The  sliding  mode  controller  (SMC)  is  a control  algorithm  that  con  offer  excelleni 
trajectory  tracking  wiihout  Ihe  need  for  precise  system  modeling.  In  many  real  world 
applicalions.  il  is  noun  easy  ia.sk  lopieci.scly  model  or  define  the  system  parameters.  The 
sliding  mode  coniroller  can  also  compensate  for  eslenai  disturbances  like  cutting  forces 
and  friction.  Tic  ihcoty  of  the  sliding  mode  controller  first  appeared  in  the  former  Soviet 
Union  in  the  1950s.  Since  then  several  researchers  have  modified  il  to  improve  its 
performance. 

Sliding  Mode  Theory 

The  single-inpul  single-output  (SISO)  sliding  mode  controller  described  here  is 
based  on  Ihc  work  of  Sloline  [17j.  Consider  the  system  in  Equation  5-1 


= f(,X)*BWu 


(5-1) 


where  1 Is  Ihe  desired  ampul,  ii  is  the  control  input,  X is  die  suile  vector.  represents 
the  input  gain  and  J{X)  lepresenus  Ihe  system  dynamics.  If  f(,X)  is  the  estimated  system 
dynamics  then  Equation  5-2  represenu  the  uncertainty  in  Ihe  system  dynamics 

^-,/l  = F (5-2) 

where  F is  the  upper  boundary  for  the  system  uncenainly. 

The  purpo.se  of  this  controller  is  to  force  the  stale  X to  track  a desired  stale  Xj  in 
the  presence  of  disturbances  and  system  uncenainlies.  Equation  5-5  is  the  tracking  error 

X = X-Xj  = [i.s (5-31 

Furthermore,  ihe  sliding  mode  controller  requires  high-speed  switching  of  the 
control  effort  to  drive  the  plant  on  a time-varying  surface  S(i)  in  state  space  R''”. 

Tracking  X=Xj  is  equivalent  to  that  of  remaining  on  the  lime  varying  surface  (or  sliding 
surface)  for  all  l >0,  Although  there  are  several  ways  to  represent  the  sliding  surface,  the 
most  common  representation  is  Equation  5-4 

+ (5-4) 

where  sf.r.r)  is  the  scalar  representation  of  the  sliding  surface.  X Is  a user  chosen  parameter 
gieaier  than  zero,  tmd  the  term  n is  the  system  order. 

For  example.  Equation  5-5  represents  a second  order  system 

i=/[j.x)  + n (5-5) 


If  n s 2 then  the  sliding  surface  is. 


5{jr.  0 = x + Xj 


(5-6) 


Figure  5-1  shows  ihe  sliding  surface  on  ihc  phnse  plane.  The  straight  line  in  the  Figure 
represents  Equation  5-6  where  \ is  the  slope  of  the  line. 

Equation  5-7  reprosenu  thedenvative  oflhe  sliding  surface 

i ‘ x-ii  + Xx  = /*u-xj*\x  (5-7) 

The  derivative  of  the  sliding  surface  is  also  called  Ihe  average  of  the  system  dynamics  on 
both  sides  of  the  surface.  When  the  system  is  in  sliding  mode  the  dynamics  can  be  written 

r = 0 (5-8) 

From  Equation  5-7  and  Equation  5-8,  the  approximated  control  signal  is 

u = -/+i,i-Xx  (5-9) 

In  order  for  the  system  to  remain  on  the  surface  the  following  Lyapunov  function 
must  be  satisfied 

i|^r^S(-ri)|r|  (5-10) 

where  7)  is  strictly  a positive  constant.  The  equation  stales  that  Ihe  squaiod  distance  to  Ihe 
surface,  mea.sured  as  s^  decreases  along  nil  system  trajectories.  This  cnsuies  that  all 
trajectories  point  towards  Ihe  surface  5(r). 

A discontinuous  term  must  be  added  to  u to  satisfy  the  Lyapunov  function. 
Equation  5-11  repiescnis  the  dlsconlinuily 


Figure  5‘L  Phase  pl^e 


i-tx  5gn(i) 


where  &gn  is  the  sigmim  runclion: 
if5>0 
sgnls)  = - I if  s < 0 

The  puramecer  k - i(x,i)  should  be  large  enough  to  guarantee  ihac  Equation  S-ID  is 
satisfied.  Ftora  Equation  5-1 1,  u isreplaced  in  Equation  5-9.  Then  the  result  is  combined 
with  Equation  5-7  to  yield  the  following  equation 

s = l/--/-*5gnU)l  (5-12) 

Subsliluiion  of  Equation  5-12  into  the  Lyapunov  function  Equation  5-10  yields 
Equation  5- 1 3. 

iwf-t-n  (5-13) 

Then  the  control  law  is  derived  in  Equation  5-14. 

II  = -/+Jj-^-i(Xsgn(j)  (5-14) 

A problem  with  the  control  law  in  Equation  5-14  is  that  it  is  discontinuous  along 
the  .surface.  The  high-speed  switching  cuuses  the  sysiem  to  chatter  along  the  surface. 

This  chattering  leads  to  undesirable  high  control  activity  and  excites  high  frequency 
dynamics  in  the  system  thut  were  neglected  in  the  model. 

A thin  boundary  layer  along  the  switching  surface  5(r)  is  used  to  smooth  the 
discontinuity.  Figure  5-2  shows  the  boundary  layer  along  the  sliding  surface.  The 
boundary  layer  is  described  by  its  width  t ond  thickness  ^ 

To  ensure  that  ihe  signal  Is  continuous  and  inside  the  boundary  layer  the  sgn  (s) 
term  in  Equation  5-14  Is  replaced  by  the  saturation  of  s/^  in  Equalion  5-1 S. 


Figure  5-2.  Boundary  layer 
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(5-15) 


saiM  = y 


if(lytsi) 


sai{y)  - sgn(y)  otherwise 

The  ihicknetts  of  the  boundary  layer  y clianges  with  time,  but  in  this  case  it  is  set 
as  a constant  parameter  to  simplify  the  calculations. 


A sliding  mode  controller  will  be  further  developed  for  the  bullscrew  drive  case. 
Figure  5-3  shows  the  open  loop  block  diagram  of  the  ballscrew  drive.  The  system  was 
modeled  for  only  a second-order  system.  The  dynamics  due  to  the  ballscrew  Hexibiliiy 
were  not  Included  in  the  model  and  are  considered  as  an  uncertainty. 

The  control  signal  ii  is  the  command  from  the  SMC,  while  r is  Ihe  position 
measured  in  counts.  is  the  disturbance  torque  coming  from  Ihe  culling  forces  and 
friction.  Kg  and  are  gains,  J is  the  inertia  of  the  screw  and  motor.  B is  the  viscous 
fricbon  in  the  motor,  and  r,  is  the  ballscrew  ratio.  Equation  5-15  represents  ihe  algebraic 
form  of  the  above  block  diagram. 


Applying  a Sliding  Mode 


(5-16) 


Rearranging  Equation  5-16  results  in  the  following  equation. 


7^(r)4fl^(r)  = u{i)-d(n 


The  acccleraiion  in  the  drive  can  then  be  expretaed  as 

'x  = j{-B,k*u-d)  15-18) 

The  exlemal  disturbance  d caused  by  the  culling  forces  and  the  friction  is  considered  a 
known  value  that  remains  between  upper  (d*)  and  lower  (d*)  limits  which  are  measured  on 
the  machine.  The  di.srurbance  can  ihen  be  tracked  using  the  following  observer  function 
d(4)  = d(t-I)  + picsr  (5-191 

where  d is  the  estimated  disturbance,  p is  an  adaptation  parameter  (p  s|0.,50).  Tis  the 
uunirol  period,  and  k is  used  to  Impose  limits  on  the  disturbance, 

0 if(d(4)Sd')  andlSSO)] 

0 if(d(4)2d*)  and[SS0)^ 

I otherwise  J 

The  sliding  surface  is  defined  as  a second  order  system  and  the  control  low  is  defined  by 
Equation  5-21 

“ = 2,lX(i,y-i)  + irfi  + B^t  + d-h*sar(i)  (5-21) 

Figure  5-4  shows  the  block  diagram  of  the  sliding  mode  controller  as  it  is  implemented  on 
the  HSMl.  A computer  simuiotion  using  SIMULLNK  shows  the  commanded  and  actual 
simulated  positions  in  Figure  S-S,  The  50,800  count  command  is  equivalent  to  a move  of 
2 inches.  The  feedrate  is  set  at  200  in/min  and  the  acceleration  isselat  I g.  Tlic  value  for 
X is  chosen  arbitrarily  while  the  value  of  4 is  calculated. 


Amplilier  Cn 

Figure  F-3.  Simplified  model  of  the  axis  dynamics 


Figure  5*4.  Block  Diagram  of  lire  Sliding  Mode  Comrullcr 


Figure  S*S.  Commanded  and  aciual  aimulaied  positions 
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The  maximum  error  beiween  (he  commooded  and  actual  simulated  poaiiions 
occurs  during  Lhe  accderaiion  period  and  is  equal  to  60coums. 

Figure  S*6  shows  the  commanded  and  aciual  simulated  velocities.  The  SMC  kept 
The  error  between  the  commanded  and  actual  simulated  velocities  to  a minimum  of 
2,000  counl/s. 

The  SMC  is  then  implemented  on  the  y-axis  of  HSMt.  The  tachometer  and 
encoder  are  used  for  feedback.  The  controller  is  used  to  generate  the  position,  velocity 
and  acceleration  commands.  Figure  5-7  shows  the  results  from  a move  of  SO.SOO  counts. 
The  system  is  stable.  A variety  of  position,  velocity,  and  acceleration  cotttmands  were 
tested,  and  they  all  resulted  in  stable  respon-ses. 

The  maximum  etror  between  the  simulated  and  actual  position  occurred  during  the 
acceleration  period  of  the  move.  Figure  5-8  is  a close-up  of  the  commanded  and  actual 
positions.  Although  the  150  count  error  during  the  deceleration  phase  is  smaller  than  that 
in  the  acceleration  phase  it  is  mors  dcinmcnial  since  it  will  cause  an  overshoot.  In  order  to 
remove  some  of  the  error  different  values  for  k and  A were  tested,  but  none  of  them 
yielded  belter  resells.  The  magnitude  of  tire  error  decreased  with  lower  acceleration 
commands. 

The  drive  was  estimaled  as  a second-order  system  to  simplify  the  calculations,  this 
simplification  causes  SMC  to  minimize  velocity  and  positton  errors  but  not  acceleration 
errors.  A third  order  system  will  minimize  lhe  acceleration  error  and  may  offer  better 
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The  bandwidth  of  Ihe  system  is  65  Hz.  which  is  lower  than  the  100  Hz  achieved  on 
this  drive  previously  using  the  PID  with  di|ilal  filters.  Increasing  the  A term  to  increase 
the  bandwidth  of  the  drive  results  in  chattering. 

Figure  5-9  shows  the  cominonded  and  actual  vclocitita  measured  on  the  drive. 

The  masimum  error  is  measured  at  7,000  count/s,  more  than  three  limes  that  of  the 
simulation.  The  acrual  velocity  measured  from  the  tachometer  also  has  more  noise  than 
anticipated. 

Figure  5-10  shows  the  difference  between  the  impiemented  SMC  and  the  existing 
PID  tuning  on  the  machine.  The  SMC  slides  towards  tite  commanded  position  minimizing 
the  error. 
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Figure  3-10.  SMC  vs  PID  loop  tuning 


CHAPTER  6 

TUN[NG  LINEAR  MOTOR  DRIVES  FOR  WIDE  BANDWIDTH 


Linear  motors  do  not  have  the  same  dynamic  limiiaiions  that  ballscmw  drives 
have.  Linear  molors  are  tiol  limited  by  lowiy  damped  second-order  natural  modes  that  are 
associated  with  the  ballscrcw  mechanisint  Iherelore.  iinear  motors  are  capabie  of  high 
speed  and  high  acceicration  operation.  PID  cuntrullers  are  usually  sufficient  for  providing 
high  bandvridihon  iinear  motors.  Since  iinear  motois  are  direct  drives,  they  produce  high 
acceieraiionst  however,  diis  also  makes  them  sensitive  to  external  disiutbancc.s.  The  two 
disturbances  that  act  on  linear  molors  in  milling  are  the  friction  coming  from  the 
guideways  and  fomes  coming  from  the  culling  tool  engaged  in  the  maicrial.  'Hie  linear 
motor  is  also  sensitive  to  changes  in  mass  and  may  become  unstable  with  substantial 
changes.  Modifications  to  ihe  PID  controller  like  friction  feedforward  and  high  gain  can 
overcome  external  disiurbances.  but  the  system  lacks  (he  robustness  to  compensate  for 
mass  changes.  With  all  these  problems,  a PID  was  still  acceptable  in  this  case  because  llie 
axis  selected  for  this  analysis  has  a fixed  mass  and  the  culling  forces  in  this  direction  are 
minimal. 

AilheMTRC,  two  brushless  DC  ironicss  linear  molors  have  been  incotpornicd  on 
ihc  Z'OXis  of  the  KSMI.  Figure  6-1  sliows  lire  Z-axis  mounted  on  two  pedestals.  Each 
pedestal  has  a linear  motor,  a lincarenctHler,  and  a guideway,  AH  the  components  arc 
mounted  nexilo  each  other  on  the  top  surface  of  each  pedestal.  The  Z-lable  which  carries 
the  AB-head  and  spindle  is  mounted  on  anti-fricljon,  high-stiffness  roller  bcanng 
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Figure  6-1.  Linear  motors  on  HSM I Z-axis 
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cuniages  that  tide  on  these  guideways.  The  motor  conflguniiion  is  such  (hot  the  magnet 
assembly  is  stretched  along  the  entire  length  of  the  travel  and  are  stationary.  The  short 
rotor  is  mounted  to  the  moving  table  and  is  constrained  to  move  between  the  magnet 
trachs.  The  encoder  is  mourned  nest  to  the  motor  and  special  care  in  the  design  was 
taken  to  eliminate  the  use  of  any  flesiblc  componems  between  them.  A more  detailed 
description  of  the  linear  motor  feed  drive  system  follows. 

Anolysis  of  linear  motor  drive 

Each  of  the  Trilogy  LM4I0-4  bm.shles.s  DCIinear  motors  provided  a peak  force  of 
2.000  N.  The  mass  of  the  table  and  spindle  was  about  900  Kg.  The  AB*hcad  was  not 
mounted  on  the  Z-asis  during  this  analysis  of  the  feed  drive.  The  maximum  possible 
acceleration  was  about  0.46  g. 

Two  Advanced  Motion  Control  (A.M.Ci  amplifiers  arc  used  to  power  the  linear 
motors.  These  PWM  amplifiers  were  sijed  to  provide  the  maximum  motor  thrust  force. 
Additionally  the  PWM  frequency  was  high  enough  iiotm  inierfcn:  with  the  motor  control. 

Two  Heidenhain  LSI 83  incremeniaJ  linear  encoders  provided  the  feedback  - one 
for  each  motor.  The  encoders  resolution  is  4 counts/',  m.  The  encoders  are  also  limited  to 
a 500  mm/s  maximum  velocity  which  Is  the  maximum  velocity  of  the  Z-axis.  Software 
compensation  tables  are  inxmllcd  on  the  Z-axis  to  compensate  for  any  yaw  between  the 
two  cncodets;  this  prevents  the  two  motors  from  fighting  each  other  and  burning  out. 

The  same  Delta  Tnu  controller  used  for  the  ballscrcw  drive  was  used  to  control  the 
Z-axis,  except  no  additional  DSP  board  wax  added.  The  controller  sample.s  the  encoder 
and  uses  the  same  signal  to  provide  the  position  loop  feedback  and  after  the  signal  is 
differeniiaied  the  velocity  loop  feedbaek.  A velocity  loop  signal  is  essential  in  uny  feed 


dnve  servo  conirol.  The  si^nui  Is  used  lo  provide  udequuie  damping  for  the  system  to 
prevent  the  drive  from  oscillating  around  the  commanded  position.  Equation  6-1  shows 


how  the  signaJ  is  generated  by  differencing  two  consecutive  position  samples  and  dividing 
them  by  the  sampling  time. 


A block  diagram  of  the  PD  iProponional-Derivative)  control  system  of  the  Z-asis 
Is  shown  In  Figure  6-2.  This  block  diagram  describes  the  diffeient  components  on  the 
feed  drive.  The  block  diagram  In  Figure  6-3  replaces  each  block  element  with  the 
appropriate  transfer  function.  Foresampic,  the  'table'  block,  which  represents  the  moving 
mo-ss  of  the  Z-axis  inble,  is  modeled  as  a simple  second-order  system  with  an  acceleration 
that  is  proportional  lo  the  force  applied  to  the  table.  This  block  diagram  has  boih  discrete 
and  continuous  block  elements.  Either  the  discrete  elements  must  be  convened  to 
continuous  elements  or  the  continuous  lo  discrete  to  be  able  to  analyze  the  system. 

Analysis  of  the  block  diagram  requires  expressing  die  discieie  time  elements  as 
their  equivalent  continuous  functions  os  shown  in  Figure  6-4.  This  reduces  the  entire 
block  diagram  to  a coniinuou.s  transfer  function  shown  os  Equation  6-2. 


The  transfer  function  indicates  thul  the  table  re.sponsc  depends  on  the  mo.ss  M of 
the  table,  the  servo-period  7^.  the  proponionai  gain  ik,,,  and  the  denvalivc  gam  Kj.  The 
servo-period  is  set  to  221  jis,  and  the  mass  is  a constant  value.  The  only  variables  are  the 
proportional  gain,  and  the  dcrivniivc  gain. 


(6-1) 


2I4A:„ 


2I4A:„ 


(6-2) 


Continuous  System 
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Discrete  Sysicin 


Fi^re6-2.  Block  diugium  of  Ihc  Z-axis  posiiion  loop 


Continuous  system 


Figure  6*3.  Discrete  time  block  diagrum  of  the  Z-axis  position  loop 


Figure  6-4.  Coniinuous  time  block  diagram  of  the  Z-axis  posiiion  loop 


79 

Tuning  (he  FD  condoller  requires  an  Inilial  guess  for  and  Kj.  Equation  6-2  can 
be  wrillcn  in  the  second'Order  form  of  Equation  6-3  where  o^  is  (he  narurol  frequency  and 
^ js(he  damping  radoof  (he  system  response.  Funhemtore,  if  is  defined  in  terms  of  the 
rise  Lime  as  shown  in  Equation  6-4.  then  K-  and  ff^can  be  expretssed  as  characteristics 
of  (he  response  to  a step  input  as  shown  in  Equation  6-J  and  Equation  6-6. 


(6.4) 


(6-5) 


Kj 


2;u), 

I 676*r,r, 


(6-6) 


With  a desired  damping  ratio  of  0.7.  and  10ms  rise  time  the  gains  are  derived  and 
applied  to  the  system.  The  initial  gains  were  too  high  and  caused  the  motors  to  emit  a 
high'pitched  hum.  Iteratively,  the  rise  time  was  reioxed  while  the  damping  mtio  remamed 
constant.  After  the  hum  had  been  attenuated  with  aSOO  Ha  low-piss  fiiterandthe 
position  loop  had  been  tuned,  the  transfer  function  of  the  system  was  measured.  From 
Figure  5-5  the  resultant  transfer  function  shows  a 70  Hr.  bundwidih. 

The  high-pitched  hum  emitted  from  the  motor  occurs  when  the  gain  in  the  velrjcity 
loop  is  too  high.  The  velocity  loop  feedback  is  being  generated  digitally  from  the  position 
loop  feedback.  This  calculation  introduce  noise  into  the  system.  Priichow  [13]  suggests 
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Figure  6-5.  Actual  (solid  line)  and  modeled  (dashed  line)  Closed  loop  frequency 
response  of  die  Z-Axis  using  linear  motors. 


two  methods  to  remove  this  noise.  The  first  suggestion  is  to  oversantple  the  feedbuck 


signoJ  end  use  averaging  filters  to  smooth  the  velocity  loop  feedback.  The  second 
suggestion  is  to  introduce  a Ferraris  feedback  sensor  for  the  velocity  loop. 

The  first  suggestion  was  evaluated  on  this  drive,  but  it  did  not  produce  favomble 
results.  Oversumpling  on  the  Delta  Tau  controller  was  achieved  by  increasing  the 
sampling  lime  for  the  position  loop  and  keeping  the  sampling  lime  for  the  velocity  loop 
fixed.  This  result  improved  the  performance  of  the  velocity  loop  but  it  draslicniiy 
deieriomicd  the  performance  of  the  position  loop. 

The  second  suggestion  requires  the  introduction  of  on  accelerometer  to  the  system. 
This  accelerometer  meusures  the  relative  velocity  between  the  moving  and  fixed  pans  of 
tile  drive.  This  accelerometer  Is  commonly  referred  to  ns  Ihe  Ferraris  sensor.  The  relative 
acceleration  measured  by  the  sensor  is  then  integrated  to  provide  a velocity  signal. 
Typically,  the  velocity  signal  has  belter  resolution  and  does  not  suffer  from  noise 
associated  with  differencing  a signal.  Using  this  sensor  on  this  linear  drive  is  not  possible 
because  of  .space  limitations. 

Currently,  the  velocity  loop  i.s  providing  enough  damping  ax  shown  in  the  step 
response  of  Figure  6-6.  The  system  exhibits  a damping  ratio  of  D.S  and  an  overshoot  of 
The  response  shows  that  the  motor  is  suffering  from  an  external  disturbance,  in  this 
case  it  Is  friction  from  guideway  covers. 

External  Disturbances  in  Linear  Motor  Drives 

Up  to  this  point,  the  linear  motor  drive  could  position  to  within  20  counts  (btim), 
Reduction  of  the  steady-state  error  was  achieved  by  the  addition  of  Integral  gain  to  the 


IT.  The  integral  gain  successfully  dropped  the  steady-slate  error  to  one  count.  In 


Figure  6-6. 


efiheliri 


addlllon,  integral  gain  provided  niueh  needed  uifTness  to  the  drive.  Slatie  siifTncs.s  is 
measured  by  applying  a constant  load  to  die  drive  while  the  actual  position  is  sampled. 

Figure  6'7  shows  the  PD  controller  under  u 200  N static  force.  The  force  was 
increased  from  0 N to  200  N over  a period  of  one  second.  The  PD  controller  shows  an 
eriorof  2.8  pm,  which  equates  to  a static  stifTncss  of  0.014  pm/N.  Under  PID  control,  the 
error  was  eliminnted  which  leads  to  the  conclusion  Lhal  ihe  static  stiffness  under  PID  is 
infinite  after  after  the  response  settled  and  as  long  as  the  force  is  within  the  thrust  capacity 
of  the  motor.  Unfortunately,  integral  gain  cun  also  lead  to  drastic  errors  during  motion. 

the  actual  position.  Il  is  most  useful  forcorreedng  errors  at  the  etid  of  each  motion. 

As  an  example,  the  linear  drive  was  commanded  a motion  under  PD  and  PID 
control.  Figure  6-8  shows  the  paching  error  which  l.s  the  difference  between  the 
commanded  and  actual  position  of  the  drive.  The  errors  meosuiud  in  the  PD  loop  are  all 
constant  and  they  ore  mainly  due  to  a lime  delay  between  the  commanded  motion  and  the 
ncrual  motion.  This  lime  delay  does  not  show  up  as  an  error  on  the  pans,  it  only  means 
lhal  the  tool  will  follow  the  commanded  path  but  arrive  al  the  end  late.  Under  PID  control, 
the  time  delay  is  not  constant  anymore.  The  integral  gain  senses  that  there  is  a huge  error 
and  attempts  to  correct  for  il.  This  correction  causes  a distortion  in  the  motion. 

During  machining  several  axes  of  the  machine  move  simultaneously.  If  each  one 
of  them  has  a different  time  delay  this  results  in  distorted  motions.  Circles  become  ovals 
and  corners  become  rounded.  A fixed  and  equal  time  delay  on  all  iixes  is  necessary  to 
produce  good  parts.  Since  Pic  muchinc  has  diflcrenl  axes  dial  arc  made  from  difftrenl 
motors  and  amplificri,  it  becomes  a challenge  to  linearize  and  match  all  lime  delays. 
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Figure  6-7.  Response  of  the  linear  drive  [o  a 200  N force  disiurbnncc 
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Figure  6-8.  Commamiul  end  uciunl  irecklng  ermn  using  bolh  PD  and  PID  comrol 


Minimizing  the  time  delay  or  removing  it  completely  is  a preferred  option.  One  method  of 
removing  the  following  error  or  lime  delay  is  feedforwarding  the  commanded  motions. 
Implemenlaiion  of  a feedforward  filter  l.s  given  in  the  following  chapter. 


CHAPTER  7 

FEEDFORWARD  FILTERS 

CNC  machining  of  complex  gcomeiriex  involves  ihesimulinneous  motion  of 
several  axes.  In  Uic  case  of  HSMl.  ii  may  involve  ihe  motion  of  ail  five  axes.  Each  of 
these  axes  has  difTercm  mechanical  and  electrical  propenies.  The  X and  Y arc  ballscrew 
driven,  the  Z is  linear  molor  driven,  and  the  AB  head  uses  rotary  motors  and  worm  gear 
transmission  for  both  axes.  Each  of  these  Five  axes  performs  differently;  but  when  they 
work  together  they,  need  to  behave  similarly  to  reduce  part  defects.  In  point-to-point 
motions,  the  tracking  error  doe,s  not  cause  any  actual  machining  errors,  but  in  multi-axis 
machining,  the  difference  in  following  errors  between  axes  causes  part  defects. 

For  example,  when  machining  a circle,  both  feed  drive  axes  are  commanded 
through  asine  wave  profile  with  one  wave  shifted  by  n/2.  At  any  given  time.  Ihepusiiion 
of  each  axis  must  be  on  the  wave  to  perform  a perfect  circle.  Figure  7-1  shows  a perfect 
circle  as  a solid-line.  When  the  axes  have  differem  delays  dicn  the  circular  becomes 
distorted  as  shown  as  u dolled  line.  To  correct  the  disionion  requires  that  the  tracking  error 
on  both  axes  be  the  same.  With  very  different  axes,  this  is  not  possible  because  they 
behave  dilTereniJy.  Feedforward  filters  are  suggested  for  removing  these  errors. 

Feedforward  fillers  have  been  the  subject  of  several  researchers.  Tomixuka  [6] 
and  Week  [S|  both  investigated  the  subject.  Their  investigations  were  discussed  in 
Chapter!,  These  filters  essentially  modify  the  command  to  counter-net  svstem  dynamics 
and  feedforward  it  to  the  system  in  order  to  reduce  the  amount  of  the  tracking  error. 


X Posilton 

Figure  7.1.  Circular  moiion  with  two  dirrcrem  time  delays  in  the  Y-diiecUon 


Keedforward  Filter  Anal>!tts 


The  purpose  of  the  feedforward  filler  isio  minimize  the  tracking  error  between  the 
commanded  moiion  and  the  ariual  motion  of  the  feed  drive.  Figure  7-2  shows  a block 
diegramof  nservomechanism  under  PD  control.  The  tracking  error  is  £(z).  Theobjeciive 
of  the  feedforward  filter  ff^z)  is  to  reduce  the  tracking  error  to  zero.  Equolion?-!  is  the 
algebruic  equivalent  of  the  block  diagram  in  Figure  7-2. 

£(z)  = XU)-Hli)l«,£(z)-£,,G(z)y(:)-i-ffj(jWj)|  = 0 (7-1) 

Solving  for  the  filter  will  result  in  Equation  7-2 

H/z)  = /r‘(z)  + ff.C(:)  (7-2) 

The  feedforward  filter  is  the  inveracof  the  open  loop  transfer  function  Hiz)  and  the 
velocity  loop  feedback.  This  assumes  that  Equation  7-3  is  the  representation  of  a 
servomechanism  drive.  Equation  7-4  is  the  discrete  time  equivalent  of  a differentiator 
using  Tustin's  method. 

HU)  = \kT-  (7-3) 

2 (J-1)' 


G(z)  = (Tustin's  method)  0 

Substituting  Equation  7-3  and  Equation  7-4  in  Equation  7-2  results  in  Equntion  7-5. 

where  Al^^and 


K^jl=  l/fKT,-) 


K,g  = K, 


Figure  7-2.  Dtscrele  time  PD  with  feedforward 
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Figure  7-3.  Feedforward  filter 
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The  Iwo  lerms  A'^i^aiid  K^jf  are  reTerred  ie  a&  the  ecceletuiion  and  vetociiy  gaint^  of 
ihc  feedforward  filler  The  accelcrttiion  gain  modiric%  Ihc  occcjcraiion  pnn  of  rhe  profile 
und  ihe  velocity  lerms  modify  the  velocity  pen  of  the  profile.  For  programming 
Equation  7-6  in  ihcconlrol  loop.  Ihe  acceleration  and  velocity  are  calculaied  from  the 
trajectory  generator  fngure  7-3  shows  how  the  two  terms  are  calculated  from  the 
command. 


Fecdfurw'ard  on  I.eadscrew  Drives 

The  closed  position  loop  of  the  ballscrew  drive  is  shown  in  Figure  7-4.  The 
feedforward  filter  Is  installed  tthead  of  Ihe  position  loop  and  alters  the  command  to 
produce  the  desired  motion.  The  values  for  the  velocity  and  acceleration  feedforward 
gains  are  calculaied  and  applied  to  the  ballscrew  drive. 

The  tracking  error  wits  reduced  on  a simulated  system.  Figure  7*5  shows  the 
tracking  error  of  a PD  system  without  Ihc  filler.  The  error  is  750counls.  With  Ihe 
feedforward  filler,  the  error  drops  to  one  count  except  during  the  ucccleraiion  phase  of  the 
motion  where  it  is  about  5 counts  as  shown  In  Figure  7-6. 

This  same  filler  was  applied  to  Ihe  ballscrew  drive  of  the  machine  and  the  results 
were  similar.  The  filter  reduced  the  error  to  less  than  one  count.  Figure  7-7  and 
Figure  7-8  show  the  errors  on  the  ballscrew  drive  of  Ihc  Y-axis  on  HSMI. 

Similarly  Ihc  feedforward  filters  are  applied  to  the  linear  motor  drive  of  HSMI 
The  filters  were  used  with  a PD  conlrollcr  and  a PID  controller.  Figure  7-9  shows  how 
the  reduction  in  tracking  error  by  using  feedforward  filters  on  the  linear  motors.  The 
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Tims{sec) 

Figure  7-6.  A more  dcmilai  look  of  ihe  itnulaled  tracking  error  with 
feedforward  fillers 
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Figure  7‘7,  Aciuul  trucking  error  on  Ihe  baJlscrew  drive 


Figure  7-8.  A cioser  iook  on  Ihe  actuai  tracking  error  of  the  haii.vcrcw  drive 


Figure  7-9.  Tracking  error  for  the  linear  motors  using  PD  ujid  PiD  conTrollers 
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irdcking  error  during  steady  velocity  is  I^m  (4  counts)  using  a PD  conirollei  and  is  a 
maximum  of2  [8  counts!  in  the  acceleration  part  of  the  motion.  With  P!D,  the 
tracking  error  is  D.S  pm  (2  counts)  during  steady  velocity  and  is  2pm  (8  counts)  during  the 
acceleration  part  of  the  motion. 

.Summary 

In  summary  the  feedforward  filler  reduces  the  tracking  error  tremendously. 
Therefore,  the  corrected  multi-axis  motions  will  not  create  machining  significant  defects 
due  to  time  delays  between  the  commanded  and  actual  motion. 


CHAPTER  8 
CONCLUSIONS 

The  seleclion  of  a lineor  motor  drive  over  a boiHcrew  drive  simply  bccousc  linear 
motors  can  achieve  higher  bandwidth  with  a RD  controller  la  not  good  design  practice. 
There  are  advantages  and  disadvantages  in  both  designs  and  proper  selection  should  be 
based  on  the  evaluation  of  all  the  litnilations  as  they  apply  for  each  particular  machine. 

On  one  hand,  ballscrew  drives  have  limited  speeds  and  accelerations,  but  with  digital 
filters  or  a more  robust  controller,  the  acceleration  can  be  increased  by  wideaing  the 
bandwidth.  On  the  other  hand  linear  motor  drives  can  achieve  high  bandwidth  with  a 
simple  PID  controller  but  suffer  from  sensitivities  associated  with  direct  drives  like 
changes  in  mas.s  and  cstemal  Tomes. 

Ballscrew  Drives 

The  major  disadvantage  in  a ballscrew  drive  is  its  limited  bandwidth  due  to  lowly 
damped  second-order  dynamics  associated  with  the  ballscrew.  These  dynamics  can  be 
allenuated  with  digital  filters  inserted  in  the  velocity  loop  of  a PD  controller.  These  filters 
are  complicated  todesign  and  implement  because  they  have  to  match  system  dynamics.  If 
for  any  leuson  the  system  dynamics  change  the  digital  filters  wilt  not  be  able  to  attenuate 
the  new  dynamics;  and  furthermore  the  filters  will  act  as  a source  of  instability. 

A more  robust  sliding  mode  controller  can  achieve  a higher  bandwidth  without 
accurate  modeling.  The  sliding  mode  controller  can  cope  with  any  unceriomtics  within  its 
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boundary  layers.  The  disadvanlagc  of  this  conirollcr  is  selling  aueplable  boundary 
layers.  These  layers  are  a compromise  between  unccrtainlies  and  pcrrormancc.  Relasing 
ihe  boundaries  for  uncertainlics  decreases  performaiKe.  Conversely,  lighlening  the 
boundary  layers  increases  performance  bui  may  cause  syslem  inslabiiily. 

Linear  Molor  Drives 

The  advantage  in  using  a linear  molor  drive  is  lhal  Ihey  are  easily  tuned  lo  achieve 
a high  bandwidth.  Linear  molars  can  reach  speeds  and  accelerations  dial  surpass 
bailscrew  drives  because  they  are  kinemaiically  simpler.  Tire  di.sudvamage  of  linear 
motors  is  Ibai  they  ore  sensilive  lo  disturbances.  The  performance  of  a FID  controlled 
linear  motor  drive  U incapable  of  coping  with  large  changes  In  the  moving  mass. 
Sufficienl  changes  may  lead  lo  Inslabiiily.  Another  disadvantage  lo  linear  motor  drives  is 
the  effects  of  cutting  forces  and  friction  because  these  are  not  suppressed  by  a 
transmission  ratio.  Because  these  forces  ate  generally  not  constant  hut  are  continuously 
applied,  Ihe  continuous  force  of  the  motor  must  be  sized  to  be  greater  than  the  sam  of  the 
culling  forces  and  friction. 

Additionally  both  linear  molor  drives  and  bailscrew  drives  have  time  delays 
between  the  commanded  and  actual  motions.  These  lime  delays  are  usually  kept  constant 
lo  minimize  errors.  This  is  not  possible  for  linear  motors  because  they  need  integral  gain 
to  increase  their  sliffness  and  to  help  them  overcome  ealemel  disturbances,  integral  gain 
causes  a noa-uniform  time  delay  Feedforward  filters  arc  added  on  all  the  ices  to  remove 


the  lime  delay. 


Kccdfor'viird  l'echnlt|UM 


Feedforward  fillers  were  used  lo  improve  tracking  of  bullscrew  drives  and  linear 
motor  drives  under  PD  and  PID  control.  These  filters  were  capable  of  decreasing  mo.siof 
the  tracking  error. 


APPENDIX  A 
MACHINE  PARAMETERS 


TiUileA-1.  LmearinolDrpar;uncicri[)Sl. 

MOTOR  MODEL 
Number  of  poles 
CoillLenglh 
Coil  weight 
Megnei  iruckweighi 
Force  conslanl 
Coil  resistance  Ph-Ph 
Continuous  force  (coil  at  75  EC  rise) 
Continuous  power 
Peak  force 
Peak  power 
Peak  current 
Continuous  current 
Heal  Dissipation  constant  (no  cool) 


liM(kg) 

lb/rt<kg/M) 

Ib/A(N/A) 


tbs  (N) 


Heat  Dissipation  consiam  (air  cool)  wuils/  •>[- 


Heal  Dissipation  constant  (liquid  cool)  watts  / '>c 


Inductance  mh 

Back  emf  V/in/sec 


LM4I0-4 


14.56(369.8) 
6.50  (2.95) 
20.0(29.8) 
23.0(102.3) 
15.3 

90  (400) 

226 

450(2000) 

6700 
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Table  A- 2.  Linear  motor  amplifier  [ 1 8] 

Amplifier  model  Units 

AC  supply  voltage  VAC 

Conlimiuus  DC  current  available  from  A 
internal  powcrsupply  A 

Peak  DC  current  available  from  A 

internal  supply  (Max.  2 Seconds)  A 

Internal  shunt  resistor  ohms 

Shunt  switch-on  voltage  V 

Bus  capacilance  uF 

Shunt  fuse  A 

AC  line  fuses  A 


S3QA4QAC 
480  3-Phase 

30  A for  Three  Phase  AC  Input 
15  A for  Single  Phase  AC  Input 
60  A for  Three  Ria.se  AC  Input 
30  A for  Single  Phase  AC  Input 
30  ® 50  W 
390 
I960 

3 ® 250  VAC 
15®  250  VAC 
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TahlcA-3.  Rotary  Moior  Punimcicrs  [I9| 
Description 

Horse  power 
Kilowatts 

Speed  at  rated  power 
Matt,  operating  speed 
Continuous  torque  (stall)  at  dO  "C 
Continuous  torque  (stall)  at  25  *C 
Continuous  line  current 
Peak  torque 
Peak  line  current 
Max.  theoretical  acceleration 
Torque  sensitivity  (stall)  ±10% 

Back  EMF(line-io-line)  ±10% 

Max.  linc-to-line  volts 

DC  resai2S  *C  (linc-to-linc)  ±10% 

Inductance  (linc-to-line)  ±30% 

Rotor  inertia 

Static  friolion 
Thermal  time  constant 
Viscou.s  damping 


Units  B-804B 

HP  18.6 

kw  13.9 


RPM  2000 

N-m  78.7 

N-m  83.4 

Amps  RMS  48.0 

N-m  230.0 

Amps  RMS  147.0 

rad/sec^  27.4(KI 

N-m/AinpRMS  1.646 

V/KRPM  99.6 

Volts  RMS  250 

Ohms  0.129 


kg-n,:  0.00840 

kg  30.6 

N-m  0.90 

Minutes  48 

N-m/KRPM  0.300 
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Tabic  A-4.  Rolory  Motor  Amplifier  [19] 


Description 

BDS4-255X 

Mam  DC  bus 

Minimum 

Maximum 

250  VDC 
360  VDC 

Unregulated  logic  bus  input 

14.5-26  VDC 
@ 0.25  Amps 
14.5-26  VDC 
@ 0.25  Amps 

Fan  AC  line  input 

95 -132  VAC 
@0.75  amps 

Output  current  (RMS/Phasel 
(4S  dccC  AMB) 

Continuous  (RMS) 

Peak  (2.0  sec.) 

(Fan  Cooled) 
55.0  Amps 
110  Amps 

Output  KVA  (@330  VDC  BUS) 
Continuous  (45  degC  AMB) 
Peak  (2.0  sec.) 

20  KVA 
40  KVA 

Internal  heat  dissipation 

445  Wnits 

PWM  switching  frequency 
Motor  current  riple  frequency  10% 
Resolverexcitation  frequency 

10.0  KHZ 

20.0  KHZ 

7.0  KHZ 

Form  factor  RMS/AVG 

< I.OI 

Speed  mgulntion  (Long  Term) 

0.075% /degC 

"Minimum  Controllable  Speed: 

Standard  l2'Bit  R/D  with  SOOO  RPM  Max.  Tracking  rate 

2.0  RPM 

12-Bit  R/D  with  14000  RPM  Max.  Trucking  Rate 

3.5  RPM 

1 4-Bit  R/D  with  2000  RPM  Max.  Tracking  Rate 

1/8  RPM 

14-Bit  R/D  with  3500  RPM  Max.  Tracking  Rate 

1/4  RPM 

16-Bit  R/D  with  500  RPM  Max,  Tracking  Rate 

1/64  RPM 

APPENDIX  a 

DIGITAL  FILTERING  TECHNIQUES 


This  is  a short  introduction  to  digital  liltcring.  it  is  not  intended  to  be  u 
comprehensive  description  on  digital  lilleritig.  A more  thorough  introduction  is  included 
inSlrum  tmd  Kirk  [20|. 

Digital  filtering  in  digital  signal  processors  is  based  on  processing  sampled  signals. 
Signals  are  sampled  at  equal  lime  intervals,  which  is  determined  by  a hardware  clock  on 
the  DSP  board  orexlernolly  fmm  a dilfercnt  processor.  The  sampling  frequency  is 
referred  to  as  fj.  The  ratio  of  Ihe  frequency  of  the  signal  over  the  sampling  frequency 
times  lit  is  the  digital  Iiequency  @. 

Individual  samples  can  then  be  obtained  by  indexing  Ihe  sample  number  n.  For 
example  for  a given  frequency/ond  sampled  frequency^  the  sine  wave  can  beevalualed 
from  Ihe  following. 


A filter  can  be  represented  by  an  equation  with  an  input  X(n)  and  an  output 
There  are  two  common  types  of  Elters  FIR  and  UR:  finite  impulse  response  and  infinite 


Equation  B-2descibes  the  FIR  Filter, 


Xfin)  = A(,X(«)-sA|X(rr  - l)  + „,  + b,X(n-t)  = (B-2) 


sin(e«) 


0, 1.2.3, ... 


(B-n 


104 


105 

where  t is  the  ordet  of  the  filler  and  ihe  b's  are  ihe  filter  coefficients.  If  a delay  operator,  c, 
is  defined  the  series  can  be  written  in  terms  of  this  delay  operator  as  follows, 

War)  = (bo  + A|f'  + V^  + ...+V''W2)  (B-3) 

The  transfer  function  of  the  filter,  H(z),  cun  be  written  as  the  output  over  the  input 

We)  = ^ = + + (B-4) 

The  roots  of  the  z polynomial  In  the  numerator  of  the  transfer  function  are  known 
as  the  'zeros'  of  the  transfer  function. 

The  infinite  impulse  response  filter  depends  both  on  the  input  and  past  output 
values  of  tlie  filter  to  form  the  outpuL  The  filtering  action  can  be  written  as  follows, 

Xf(n)  + ...*a„Xf(n-m)  = b„X{n)*b,X(n- \)  * ...  + b^Xin- k)  (B-5) 
where  the  terms  areas  previously  mentioned  and  the  a’s  represents  the  coefficients  placed 
on  the  past  output  of  the  filter  This  can  be  more  compactly  written  as  follows. 


There  is  a possibility  that  the  filter  can  be  unstable.  The  transfer  function  of  the 
filter  can  be  written  in  terms  of  the  delay  operator  z as  follows. 

(1  *0,2  +0,2  +...  + a„2  ) 

The  roots  of  the  denominator  polynomial  an 


ic  'poles'  of  the  filter.  The  poles 


must  lie  within  the  unit  circle  in  the  2 plane  in  order  for  the  filter  to  be  stable. 
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